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Localization properties of particles in one-dimensional incommensurate lattices without interac- 
tion are investigated with models beyond the tight-binding Aubry-Andre (AA) model. Based on 
a tight-binding t± — ti model with finite next-nearest-neighbor hopping £2, we find the localization 
properties qualitatively different from those of the A A model, signaled by the appearance of mobility 
edges. We then further go beyond the tight-binding assumption and directly study the system based 
on the more fundamental single-particle Schrodinger equation. With this approach, we also observe 
the presence of mobility edges and localization properties dependent on incommensuration. 
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The physics of quantum transport in random disor- 
dered potentials has been a subject of substantial inter- 
est for condensed matter physicists for decades. The ex- 
tended Bloch waves in a periodic lattice could undergo 
a quantum interference induced transition into localized 
states due to random disorder by a mechanism commonly 
referred to as Anderson localization pQ. Matter waves 
can also be localized in deterministic potentials that ex- 
hibit some similarities to random disorder [2j |3j [4j IS]- 
Quasi-periodic potentials, such as incommensurate lat- 
tices (the superposition of two or more lattices with in- 
commensurate periods), are notable examples and have 
been extensively studied with the Aubry-Andre model 
[2] . Such potentials have been shown to exhibit interest- 
ing quantum transport phenomena in themselves. Incom- 
mensurate potentials, for example, are theorized to have 
fractal spectrums [6]. However, it remains challenging to 
study these phenomena in solid state experiments, as it 
is difficult to systematically control the disorder in solid 
state systems. In contrast to the solid state systems, ul- 
tracold atoms loaded in optical lattices offer remarkable 
controllability over the system parameters, making it an 
attractive platform for the study of the localization of 
matter waves. Recently, Anderson localization of non- 
interacting Bose-Einstein condensates (BEC) has been 
observed in a one-dimensional matter waveguide with a 
random potential introduced with laser speckles [7 j . Sim- 
ilar experiments have also been done in quasi-periodic 
optical lattices [8j [9]. 

Localization of noninteracting particles in one dimen- 
sional incommensurate lattices is often studied with the 
Aubry-Andre model (AA) with nearest neighbor (nn) 
hopping, where one of the lattices is assumed to be rela- 
tively weak and can be treated as a perturbation . Within 
the framework of the A A model, there is a duality point, 
at which a sharp transition from all eigenstates being 
extended to all being localized occurs. However, in ul- 
tracold atom experiments, one can tune the depth of each 



lattice in a controllable way and bring the system out of 
the tight-binding regime. To explore the physics of local- 
ization for shallow lattices, it is of interest to go beyond 
the AA model and the tight-binding assumption [TP] . 

In this work, we first study the tight-binding t\ — 
model, which extends the AA model by including the 
next-nearest neighbor (nnn) hopping. The inclusion of 
the nnn hopping destroys the self-duality possessed by 
the AA model and the localization properties of the sys- 
tem become more complex through the emergence of mo- 
bility edges. We then examine the system directly with 
the single particle Schrodinger equation. We discretize 
the equation and solve it numerically without any fur- 
ther assumption. Within this formalism, we also find the 
existence of mobility edges, consistent with the t\ — £2 
model results, and we find localization properties with 
non-trival dependence on incommensuration. 

Consider diffuse, noninteracting, ultra-cold atoms in 
a one-dimensional incommensurate lattice, where the 
atoms can only move along the x-axis. The lattice po- 
tential is given by 

V V\ 
V(x) = — cos(2/clx) + —cos(2akLX + 5), (1) 

where Vq and V\ describe the depth of the primary and 
secondary lattices respectively, Ul is the wave- vector of 
the primary lattice along the x-axis, a is an irrational 
number characterizing the degree of incommensurability 
between the periods of the two lattices, and S is an ar- 
bitrary phase (in our calculations it is chosen to be zero 
for convenience, without loss of generality). When the 
depth of the primary lattice is sufficiently large as com- 
pared with the recoil energy E r = (hki) 2 /2m as well as 
the depth of the secondary lattice Vi, the physical prop- 
erties of the system can be studied with the well-known 
single-band tight-binding Aubry-Andre model: 

t(u n -i + u n+1 ) + V n u n = Eu n , (2) 
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FIG. 1: Inverse participation ratio of all eigenstates for t\ — £2 
model with a = (y/E-l)/2. The size of the system is chosen 
to have 1000 sites. The four panels correspond to £2 = 0, 
0.01, 0.05, and 0.1 respectively. (£1 is the unit for energy.) 
Darker shading corresponds to more extended states while 
lighter shading corresponds to more localized states. 



in which only the coupling between nearest-neighbors 
(nn) is retained and the incommensurate modulating po- 
tential V n = V cos(27ran). The duality point is given 
by V/t = 2. The nn hoping term, £, is determined by 
the primary potential and can be approximated by the 
expression 
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lattice potential and its magnitude can be roughly esti- 
mated by applying Gaussian approximation for the Wan- 
nier states: 



V ~ -— exp . 



(4) 



We note that V depends on Vi, a, and Vo/E r . As a naive 
extension to the AA model, we ask what will happen if 
the coupling between next-nearest-neighbors is included. 
To answer this question, we consider the model: 



^ td(Un-d + U n+d ) + V n U n Eu n 



(5) 



d=l,2 



where V n = Vcos(27rn). We solve the equation by di- 
rect diagonalization. To quantify the localization of the 
wave function, we compute the inverse participation ratio 



FIG. 2: Inverse participation ratio on the £2 — V plane for 
a = (a/5 — l)/2 based on the £1 — £2 model. The four pan- 
els correspond to four eigenstates labeled by i, with ascend- 
ing eigenenergies. Blue regions correspond to more extended 
states and red regions correspond to more localized states. 



(IPR): 
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where the superscript i denote the z-th eigenstate (or- 
dered according to energy from low to high). For spa- 
tially extended states, IPR approaches zero whereas it is 
finite for localized states [TT] . 

Fig. [l] shows the IPR values of all eigenstates as a 
function of the effective strength V of the secondary 
lattice based on the tight-binding £1 — £2 model with 
a = (V5 - l)/2 for various values of £2 (£1 is chosen 
to be unit of energy). The calculation for Fig. [l] is done 
for a system with 1000 sites in the primary lattice. For 
small values of £2 (e.g. £2 = 0.01), the localization prop- 
erties of the system have essentially the same features as 
those determined by the AA model. However, when £2 = 
0.05 or higher, AA duality is clearly destroyed and local- 
ization transitions appear to be energy dependent. For 
lower energies, the transition can appear for V < 2t\ and 
for higher energies, the transition can appear for V > 2t\. 

In order to demonstrate the dependence of the local- 
ization transition on £2 , we show the distribution of IPR 
on the £2 — V plane for four different eigenfunctions with 
a = {y/h — l)/2 in Fig. [2] For the calculation, the size of 
the system is chosen to be 40,000. At £2 = 0, the £1 — £2 
model reduces to the AA model, and from Fig. [2j one 
can see the sharp transition when V is increased across 
the duality point V — 2. However, the localization prop- 
erty of the system is greatly complicated when £ 2 is finite. 
Besides the appearance of mobility edges, our results also 
reveal that the dependence of the localization property 
on £2 is not monotonic, e.g. at fixed V < 2 when £2 is 
increased the ground state could be tuned from extended 
to localized, but further increasing of £2 could bring the 
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ground state into an extended state again. 

We infer from the results presented in Figs, [I] and [2] 
that 1) the A A duality is destroyed by having ti ^ 0; 
2) instead of the V = 2t\ dual point, the system has 
energy dependent mobility edges for ^ 0; 3) the pre- 
cise localization condition deviates up or down from the 
V = 2ti AA condition depending on the energy of the 
eigenstate and the value of t^. As illustrated by Figs. 
[I] and [2j the t\ — model itself could be of interest. 
However, for the study of localization properties in ID 
incommensurate lattices, its validity must be dealt with 
caution, especially when is not sufficiently small as 
compared with t\. The tight binding nn and nnn hoping 
integrals t± and ti can be estimated with the Wannier 
basis, which is fully determined by the primary lattice. 
One can easily estimate that when Vb = 3E r , the ra- 
tio of tijt\ is on the order of 10%. To get higher tijt\ 
ratio, one will need to tune the lattice potential shal- 
lower and should expect the tight-binding approximation 
to break down at some point. Alternatively, to study the 
interesting physics of localization in this regime, we nu- 
merically solve the single-particle Schrodinger equation 
without any tight-binding approximation: 



(- 



2m dx 2 



V(x))ip(x) = Eip(x). 



(7) 



To achieve this goal, we discretize the Schrodinger equa- 
tion in the position basis with a finite system size of 
length L = A/a, where a is the lattice constant of 
the primary lattice associated with Vb. The continuous 
Schrodinger equation is now cast into the following form: 



( O ) 



(Vq cos(2kLn5) 



2m J S 2 



(8) 



where 5 = Na/M is the step interval for the discretiza- 
tion with M denoting the total number of steps. Then we 
proceed by diagonalizing the M x M matrix of the dis- 
cretized Hamiltonian and study the first N eigenstates 
with smallest energy eigenvalues. These states would 
correspond to the ground band for the case with no sec- 
ondary lattice (i.e. V\ =0). In our calculations for the 
following results, we have set TV = 500, M = 80, 000, and 
2k L = 1. 

IPR values (obtained with Eq.([6| by replacing u n with 
ifj n ) of the first N eigenstates as a function of the sec- 
ondary lattice strength V\ are shown in Fig. J5] for a 
primary lattice strength of Vb — 30E r . In Fig. (3fa) the 
irrational ratio a is set to be the inverse golden mean, 
(y/5 - l)/2 whereas in Fig. ^b), a = tt/2. The bold- 
dashed line represents the AA duality point calculated 
with Eqs. Q and Q. We can see that the localization 
properties shown in Fig. [3] closely resemble the well- 
known results from the AA model (see top panel in Fig. 
[T]). We do note, however, that the IPR results of Fig. 
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FIG. 3: Inverse participation ratio obtained by solving 
the Schrodinger equation and calculated AA duality point 
(dashed line) at V = 30E r (a) a = (v 7 ^ - l)/2; (b) a = tt/2. 
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FIG. 4: (a) Inverse participation ratio obtained by solving 
the Schrodinger equation and calculated AA duality point 
(dashed line) at Vb = 2E r (a) a = (y/E — l)/2, solid lines are 
estimated location of localization transitions; (b) a — tt/2. 



[3] indicate a dependence on the specific value of a with 
a = (y/E — l)/2 providing a sharper A A duality than 
a tt/2. 

In Fig. |4ja) we show the IPR values for the case of 
Vb = 2E r and a = (y/E — l)/2. In this case, the eigen- 
states no longer appear to localize all at once, but in dis- 
crete steps (represented by the solid lines in the figure). 
This localization behavior is similar to what we observed 
in the t\—t2 model (see bottom panel in Fig. [TJ. Also the 
transitions occur at fairly large values for Vi, where the 
secondary lattice can no longer be treated as a pertur- 
bation. We have also studied the cases where Vb = 2E r , 
a = tt/2 (Fig. Qb)) and a = (y/E + l)/2 (not shown in 
the figure). In these cases no localization was observed 
in the eigenfunctions for any value of V\ investigated (up 
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FIG. 5: Inverse participation ratios of the ground state wave- 
function for the case Vb = V\ and a equal to fractional multi- 
ples of (y/E — l)/2. The solid blue line represents an approx- 
imate analytical boundary between localized and extended 
regions based on the A A duality point. 



to Vi = Vb). This suggests that incommensurability be- 
tween the lattices is not a sufficient condition to observe 
localization for shallow cases. 

To examine the dependence of the localization transi- 
tions on a, we set Vb = V\ and calculate the IPR of the 
ground state for various values of Vb and a (the values of 
a examined are all proportional to (y/E — l)/2). These 
results are shown in Fig. [5) We see fairly distinct re- 
gions of localized and extended states, with localization 
tending towards areas of larger values for Vb and smaller 
magnitudes for a. The blue line in Fig. [5] represents 
the set of points (a,Vb) such that the A A duality point 
(calculated from Eqs. ([3| and Q) is equal to the lat- 
tice strength Vb- These sets of points serve as a simple 
heuristic estimation of the boundary between localized 
and extended states based on AA duality condition. Al- 
though in principle we should not expect the AA duality 
point obtained from Eqs. (|3| and Q to be applicable 
in the case of shallow lattices, this simple analytical re- 
sult is in good qualitative agreement with our numerical 
findings. 

We now briefly discuss how some of these results may 
be observed in cold atom experiments. We consider a 
diffuse BEC that is loaded into an incommensurate op- 
tical lattice, confined by a harmonic trap, Vtrap = ^lx 2 . 
We assume that the diffuse gas is prepared in the ground 
state. At time T = 0, the harmonic trap is suddenly 
turned off and the BEC is allowed to diffuse. Localiza- 
tion can be observed by monitoring the IPR of the den- 
sity wave function over time. In Fig. |6j we present the 
calculated values for the IPR as a function of V\ for the 
wave function after a fixed period of time, Tq « h/E r , 
has passed since the trap was turned off for the cases with 
Vb = 2E r , ft/E r « 1(T 7 , a = (y/E - l)/2 and a = tt/2. 
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FIG. 6: Inverse participation ratio of ground state wavef unc- 
tion at time To ~ h/E r after the trap potential Vt rap = Qx 2 
has been turned off (Q/E r « 10 -7 ). 



In the figure, we see the two cases are similarly delocal- 
ized for small values of V\ . But for larger values of V\ , the 
IPR for thea = (y/E-l)/2 case begins to grow, showing 
increasing degree of localization, while in the a = ir/2 
case it remains constant. 

In conclusion, we have studied the localization proper- 
ties of noninteracting particles in a one-dimensional in- 
commensurate optical lattice system based on a tight- 
binding t\ — ti model with nearest-neighbor as well as 
next-nearest-neighbor hopping. We reveal the emergence 
of mobility edges when the next-nearest-neighbor hop- 
ping is finite. We have also gone beyond the tight-binding 
approximation by directly modeling the system with the 
fundamental single-particle Schrodinger equation, which 
is expected to provide more reliable theoretical descrip- 
tion of the system especially for the case with shallow pri- 
mary lattice potential. By diagonalizing the discretized 
Hamiltonian, we numerically solve the Schrodinger equa- 
tion. Our results clearly show the existence of mobility 
edges. Our study also reveals that the emergence of lo- 
calization is sensitive to the magnitude of the irrational 
ratio a of the incommensurate lattice potentials when the 
system is well outside of the tight-binding regime. Our 
results also establish the fragile nature of the AA duality 
which gives way to mobility edges as soon as longer range 
hopping, even at the nnn level, is turned on. It will be 
interesting to verify our predictions about the sensitive 
qualitative dependence of ID incommensurate localiza- 
tion on Vb, Vl, E ri and a through experiments in cold 
atomic systems [7J [H [9] . 
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